Science (UAEMS) 
Peer-Reviewed Journal 


INFQGAIN 


PUBLICATION 


International Journal of Advanced Engineering, Management and 


ISSN: 2454-1311 | Vol-9, Issue-2; Feb, 2023 


Journal Home Page: https://ijaems.com/ 
Article DOI: https://dx.doi.org/10.22161/ijaems.92.2 


An Efficient & Less Complex Solution to Mitigate 
Impulsive Noise in Multi-Channel Feed-Forward ANC 
System with Online Secondary Path Modeling 
(OSPM) 


Fahad Ihsan Khan!, Muhammad Umer Javed’, Didem Kıvanç Türeli’ 


!3Istanbul Okan University, Türkiye 
2NUST, Islamabad, Pakistan 


Received: 11 Jan 2023; Received in revised form: 08 Feb 2023; Accepted: 16 Feb 2023; Available online: 23 Feb 2023 


Abstract— This paper deals with impulsive noise (IN) in multichannel (MC) Active Noise Control (ANC) 
Systems with Online Secondary Path Modelling (OSPM) employing adaptive algorithms for the first time. It 


compares performance of various existing techniques belonging to varied computational complexity range 
and proposes four new methods, namely: FxRLS-VSSLMS, VSSLMS-VSSLMS, FxLMAT-VSSLMS and NSS 
MFxLMAT-VSSLMS to deal with modest to very high impulsive noise (IN). Simulation results show that these 
proposed methods demonstrated improved performance in terms of fast convergence speed, lowest steady 


state error, robustness and stability under impulsive environment in addition to modelling accuracy for 


stationary as well as non-stationary environment besides reducing computational complexity many folds. 


Keywords— impulsive noise, multichannel, Online Secondary Path Modelling 


I. INTRODUCTION 


Periodic noise, typically the low frequency noise, is a 
serious issue in many noise handling applications, such as 
those in the industry, production plants, air conditioning 
units, within aerial vessels, ships, and other vehicles [1]. It 
also imperils human mental and physical health, particularly 
for the infants and the older ones [2, 3]. Active-noise- 
control (ANC) is the most useful instrument for reducing 
this low frequency noise [4]. ANC system works on the 
superposition principle, which reduces the impact of 
undesired noise by causing destructive interference between 
the acoustic waves from the noise source and anti-noise 
signal generated by the noise controller [5]. It uses adaptive 
digital filters to track the noise source, acoustic atmosphere, 
and unknown acoustic paths [6]. 

An ANC system can be configured in feedforward 
(FF), feed backward or, combination of both, as hybrid 
design [6]. A single channel (SC) feedforward (FF) ANC 
system encompasses a reference and an error microphone, 
in addition to a loudspeaker known as secondary source, as 
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depicted in Fig. 1. The reference and error microphones are 
used to pick-up undesired x,. (n) and residual e(n) noises, 
respectively. With the help of these two signals, an anti- 
noise signal y(n) with the same amplitude as of undesired 
noise signal but opposite in phase is generated by the noise 
control adaptive filter called control filter represented by 
W(n) and released through secondary source 
(loudspeaker). The cancelling signal y(n) propagates 
through secondary path S (n) to the error microphone where 
residual error e(n) is calculated which is iteratively used to 
update control filter’s W(n) coefficients [7]. 
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Fig. 1: Single-channel (SC) feedforward (FF) ANC system 


In an ANC system, coefficients of control filters 
are updated recursively on each iteration using adaptive 
algorithm [4, 5]. Least Mean Square (LMS) algorithm and 
its different variants have proven most effective adaptive 
algorithms in ANC systems due to its simplicity and 
performance [8]. Filtered-x-Least Mean Square (FxLMS) is 
a renowned, simplest and commonly used adaptation 
algorithm for ANC systems [8, 9]. The input reference 
signal x,(n) in the FXLMS algorithm is passed through a 
prototype of the perceived secondary acoustic path $(n), 
succeeding adaptive noise controller W (n), and therefore 
termed as "filtered x algorithm". Although, FxLMS 
algorithm is moderately resilient to inaccuracies between 
secondary path S(z) and modeling filter $(z);however, the 
noise reduction capabilities are lower than to what are under 
ideal environment [9]. Thus, to cater for varying nature of 
secondary path and its effect on the overall performance of 
an ANC system, online identification or modeling of 
secondary path is essential to preserve stability, and 
robustness and to keep noise reduction performance at 
optimum level [10]. 


NOISE 
SOURCE 
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Eriksson and Allie [11] proposed a method to 
approximate the secondary path coefficients by infusing an 
internally generated auxiliary random white noise vu, (n) 
into the system as shown in Fig. 2. Modeled secondary path 
S(z) is estimated by measuring response of original 
secondary path S(z) to this auxiliary white noise vu, (n). 
Fig. 2 shows the use of an extra online secondary path 
modelling (OSPM) adaptive filter $(z), based on LMS 
algorithm, in addition to the FxLMS algorithm-based 
control filter W(z). The reference noise signal x,(m) 
propagates through primary path P(z), (between input noise 
source and error microphone) while cancelling signal 
propagates through secondary path S(z), (from secondary 
source (loudspeaker) towards error microphone) where 
residual error e(n) is calculated to update control filter 
W (z). A glaring issue that Eriksson’s method faces is that 
auxiliary random white noise v, (n) turns up in the residual 
error signal e(n) which interferes with convergence of 
control filter W (z). Similarly, the residual error e(n) itself 
causes hindrance to the modelling process and convergence 
of OSPM filter. Various researchers proposed different 
methods for improvement in Eriksson's technique. For 
example, Akhtar et al. proposed improvement in Eriksson's 
method by using a simple VSSLMS algorithm for 
modelling process [12]. Ahmed et al. [13] proposed 
betterment through two stage Auxiliary Noise Power 
scheduling (ANPS) to calculate time varying gain G(n).. 
Yang et al. [14] innovated use of self-tuning power 
scheduling for determining the gain G(n) based on the 
variations of the modeling error g (n) whereas three 
adaptive filters are used to tune the step size |i, of OSPM 


filter. Additionally, to speed up the process of reducing 
disturbances in the modelling process, an additional 
reference signal cancelling filter H(z) was also utilized in 
[14]. 


Fig. 2: Online Secondary Path Modeling (OSPM) by Eriksson 
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A single-channel ANC system (as shown in Fig. 1 
& Fig. 2) is an effective tool to manage low frequency noise 
in a narrow duct. However, in more practical scenarios, as 
the noise field grows in an extended channel or in any large 
enclosure, the process of noise control or cancellation 
becomes further intricate than in a simple narrow channel. 
Therefore, in order to manage such comprehensive and 
complicated noise scenarios, it becomes imperative to 
deploy a multichannel ANC (MCANC) system made up of 
several secondary sources, numerous error sensors, and 
different reference mics [15]. In Fig. 3, a (Ix/ x K) 
MCANC structure with I number of reference noise inputs 
xy, (n), fori=1,2,3,..... ,I, J number of secondary 
loudspeakers, and K number of error microphones e; (n), 
fork =1,2,3,..... , K, is shown. In such a construction of 
MCANC system, a total of J x K time varying secondary 
paths 5,;(n) are required to be modelled through OSPM 
process which is comparatively a more cumbersome job 
than for a single channel ANC system. This increase in 
complexity is mainly due to generation of complicated error 
signal e, (n) from each of k error microphones which is a 
combination of signals arriving from different primary 
Py;(n) and secondary S,;(n) paths. S.M. Kuo and D.R. 
Morgan concluded in [4] that different reference signals 
x, (m) may be averaged over to make single reference 
signal x,(n) which will have same effect as each of the 
individual reference signals combined but with reduced 
complexity. Thus, a general 1 x J x K MCANC system will 
be considered in this paper. 


Enclosure 


NUT 


NOISE 
SOURCE 


Fig. 3: Structure of I X J] x K MC ANC system 


Most researchers, including [12, 13, 14, 16], 
primarily dealt with Gaussian noise, however, many 
practical applications in industry and construction like 
stamping and cutting machinery in the production set-ups, 
pile drivers and gasoline engines etc. require dealing with 
impulsive noise (IN) [17, 18]. IN is often characterized by 
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strong, abrupt mutations that change the signal's distribution 
into a non-Gaussian one [19]. A symmetric a-stable 
distribution can be used to model impulsive noise with the 
following characteristic function [17]: 


face (1) 
Where a in above equation is the characteristic exponent 
that determines form of a distribution and its value ranges 
from 0 < æ < 2. Tail will be heavier for æ closer to 0, 
depicting very high impulsive nature of the noise. On the 
other end, as œ approaches 2, the impulsive nature keeps on 
diminishing and distribution becomes more and more 
Gaussian. y, in (1) is the scaling parameter which if set to 1, 
makes the SaS distribution a standard distribution. Same 
will be utilized in this paper i.e., 0 < æ < 2andy = 1. Fig. 
4 shows a standard SaS distribution for different values of 
a. 


Symmetric a-stable densities - G=0, ^21, d6=0 


f(x) 


-5 0 
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Fig. 4: PDFs for standard SaS distribution of various a 


The traditional FxLMS algorithm, which minimizes error 
signal's variance, has shown to be quite effective at 
suppressing Gaussian noise [8]. However, due to non- 
existence of second-order moment in impulsive noise, the 
FxLMS method degrades and diverges when dealing with 
IN [20]. Filtered-x least mean p-power (FXLMP) approach 
as described by Leahy [21] works on minimizing the lower 
order fractional moment of error, p (0«p«a), which occurs 
for stable distributions. In terms of active IN suppression, it 
is more resilient than the FxLMS algorithm. Its convergence 
speed, however, is very slow, particularly when noise 
becomes increasingly impulsive. Furthermore, it 
necessitates a difficult process of prior approximation of p 
based on a for adequate results. Sun [22] and Akhtar [23] 
proposed modifications in FxLMS by applying threshold in 
the algorithm. This threshold eliminate the aberrant values 
in reference noise and / or error signal/s, which may lead to 
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instability in the ANC system. Sun used Modified 
Reference FXLMS (MRFxLMS) algorithm [22] to discard 
any value of reference signal above a certain threshold while 
Akhtar presented Threshold based FXLMS (Th-FxLMS) 
algorithm [23] to clip samples of reference noise and / or 
error signal/s above threshold to provide stability in the 
system. Similarly, Zeb [24], in 2017, employed same 
threshold concept using Filtered-x Recursive Least Squares 
(FxRLS) algorithm and suggested Threshold based FxRLS 
(TFxRLS) algorithm for IN in ANC systems. The TFxRLS 
algorithm improved the convergence of the system greatly 
but conceded on increased computational complexity. Zeb 
[24] also proposed hybrid Modified Gain FxRLS- 
Normalized Step Size FxLMS (MGFxRLS - NSSFxLMS) 
algorithm that shows faster convergence speed than 
NSSFxLMS with lesser complexity than FxRLS. 
Researchers, in [22, 23, 24], worked with offline modelling 
in their threshold-based algorithms as computation of 
optimal threshold was not possible in the case of ANC 
systems with OSPM. Jabeen et al. [16] employed OSPM for 
IN in single channel ANC system using FxRLS with the 
help of an additional filter as in Yang’s method [14] and 
proposed three different variants namely, FxRLS-FxRLS, 
MGFxRLS-FxRLS and VSSFxNLMS-FxRLS. These 
variants improved modeling accuracy, provided faster 
convergence and robustness with varying computational 
complexities. Although, FxRLS and its variants [16, 24] 
improved performance, but its complexity is of order O(n’) 
as compared to the order O(n) complexity of the algorithms 
belonging to FxLMS family. 


Another technique to deal with IN, but still keeping 
the complexity low, is to use high-order-error-power 
(HOEP) adaptive algorithms [25, 26, 27]. Filtered-x Least 
Mean Absolute Third (FXLMAT) algorithm is one such 
HOEP algorithm that uses mean absolute third power of 
error signal e(n) and outperforms FxLMS for majority of 
noise probability densities [26, 27, 28]. However, FXLMAT 
faces stability issues because of its dependence on various 
input characteristics i.e., weight initialization and variance 
etc. [28]. Xiong, in [29], formulated normalized LMAT 
(NLMAT) to cater for these instability issues which can 
subdue non-Gaussian noise better than other algorithms. 
Moroever, H. A. Khan et al. [30], employed FxLMAT to 
mitigate IN in single channel ANC System with OSPM. 
Khan proposed different variants, FXLMAT, MFxLMAT 
and VSSFxRNLMAT [30] to be used in control filter W (n) 
while employing VSSLMS in modeling filter $(n) in 
addition to a third filter H (n) as in yang's method [14]. 


Various researchers [22, 24, 20, 23] presented 
different methods to alleviate IN but almost all of them 
worked with offline modeling of secondary path which are 
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not reliable solutions for time varying paths. Recently, 
Jabeen et al. [16] and Hashir et al. [30] worked on Single 
channel ANC system with OSPM for IN. However, there is 
still no published research that can control the IN actively 
and adaptively in a MC ANC systems along with the 
employment of OSPM technique. Considering this fact, we 
took motivation to undertake this challenging task and 
present out proposed methods for such a scenarios in this 
paper. 

The layout of remaining paper is that Section-II 
briefly discuss Akhtar’s method for MCANC system with 
OSPM. Section-III presents new proposed techniques, 
followed by section-IV which comprises of complexity 
comparisons of various methods discussed in this paper. 
Section-V shows computer simulations to validate results of 
proposed methods, summed up by concluding note. 


II. BASIC MC ANC SYSTEM WITH OSPM 


Akhtar et al. [12] presented an efficient VSS-LMS 
algorithm based MC ANC system for Gaussian input using 
Eriksson's [11] structure of OSPM to keep the complexity 
low (Fig. 5). This variable step size (VSS) strategy used 
power ratio p, (n) of error signals e,(n) and f(n) to 
compute step size LL, (n) for OSPM filters as follows: 


us, 0) = pxGOus min + (1 Pk) )Usmax 2) 


Where Hsmin and p,,,,, are lower and upper step size 
values, determined experimentally. Power ratio py (n) is 
given by: 

Py, (n) 
Pe, (n) 
These powers P,(n) are estimated through low pass 
estimator (LPE), as: - 


Pn) 2AB(n-1)-(1-2y?(n ; yf,e (4 


Whereas error signals are defined as: - 


px(n) = (3) 


ex (n) = dy (N) - [Y a * y',,60] 7 [va * v4] ; d 


and 

fe 0) = ex (n) — [95 00) + 900] (6) 
Akhtar et al. initialized OSPM filters by offline modeling 
which was stopped when error was lowered to -5 dB. 
Afterwards, OSPM filters are updated using following: - 


Sj; (n 1) = S, (n) + Us, M [fi G) v;(n)| (7) 


After convergence of OSPM filters, Control filters in 
Akhtar's method are updated by MeFxLMS algorithm as: - 


wj (n+ 1) = wi (0) + uv; [f G0 Bam) + fm) £500] (8) 
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Fig. 5: Akhtar's 


III. PROPOSED METHODS 
A. Proposed Method - 1: FxRLS & VSSLMS. 
FxRLS algorithm has been proven favorable in dealing with 
impulsive noise in single channel ANC system [16, 24]. In 
this method (Fig. 6), FxRLS has been employed in control 
filters of MC ANC system while modelling filters are 
adapted through VSSLMS algorithm as follows: - 


$ (n 1) = Sj Q) + us Q) [fi(n) vjQ)] (9) 
where f; (n) is modelling error signal, defined as: - 
fir) = ey (n) — [P'ra ™) 8,00] (10) 


and 
ex (n) = dy (n) - [y'a 6 + y',,00] - 
d, (n) = py, * x, (n) 


while step size 4s ,is calculated through: - 


As, (1) = Pr) Us min 2 (1 (12) 


Usmin 204 Usmay are lower and upper step size values 


determined experimentally and p(n) is power ratio of 
modelling and residual errors, given as: - 


— Pk (n) )us max 
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d4(n) 


$2) | 


method for OSPM in MC 1x2x2 ANC System 


P fk (n) 
Pe, (n) 


where Py (n) and P,, (n) are powers of f, (n) and e, (n) 


px (Qn) = (13) 


error signals which are estimated employing low pass 
estimator of the form: 


EQ) 2 AB, - 1) c (17 2v, ; y 2 f,e (04) 


After convergence of OSPM filters, FxRLS algorithm is 
used to update ANC filters as following: - 


w (+1) = w Q) + fx) Ky Q0) + fo) Kw Q0) (15) 
where 
Pyg kj) 


x Tan) Pyyj x P216) + a 


Ky, (n) E 


Pwy; is initialized as Py, (0) = 6711 and 6 is regularization 


parameter having experimentally determined value. 
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Fig. 6: Proposed method — 1, FxRLS & VSSLMS 


B. Proposed Method — 2: VSSLMS & VSSLMS. ŝkj (N+ 1) = Se; (n) + Us, (n) [f G0 vj(n)] (17) 
Proposed method - 1, FxRLS & VSSLMS, produced good Whereas control filter is updated using: - 
results but at the cost of high computational complexity. w; (n+1) = w; (n) + bw (n) [f. (n) R'is (n) + fo (n) LA (n)](18) 
Consequently, there was a dire need to find a less complex 
solution to handle impulsive noise. In this method, a less Step size Hw j (n) to be computed as under: - 
complex algorithm, VSSLMS, is used in both modeling and a= j 

= _ c (1 — pn 19 
control filters that yet shows performance comparable to Hy Q0) = PRODI i ( Prl Jonas (9) 
that of FxRLS variant discussed above. Modeling process is And py (n) is calculated using (13) and (14) 
same as given in Proposed method 1 (equations 9-14). Thus, 
following equations are used to update OSPM filters: - 
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Fig. 7: Proposed method — 2, VSSLMS & VSSLMS 


C. Proposed Method — 3: FXLMAT & VSSLMS. 
Proposed Method - 2 provided simplicity, but convergence 


wj(n +1) = wn) + n, QD [fi^ sign fA) Q0 + fr" (n)sign 


h is determined imentally. 
speed was affected. To find fast convergence along with less M prd a » e ane oes s SSLMS 
elling filters are ed throu 
complexity, a HOEP algorithm, FXLMAT in combination " ith i , : i : ie 14 : 
algorithm using equations 9- 
with VSSLMS (Fig. 8) is tried as under: 50 Seg 


Control filters will be adapted using: - 
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d4(n) 


f2(n) 
P2(n) 


7 f,(n) 
p,(n) 


Fig. 8: Proposed method — 3, FXLMAT & VSSLMS 


D. Proposed Method — 4: NSS MExLMAT & 
VSSLMS. Bona fide FxLMAT did not produced 
anticipated results. A modified version of FxLMAT is 
proposed here which has much less computational intricacy 
but shows performance comparable to that of FxRLS 
variant discussed above. 


Here, again, the modelling process is carried out using 
VSSLMS technique similar as given in Proposed method 1 
(equations 9-14). Therefore, details of only control filters 
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will be discussed here. Below is the weight update equation 
for control filters: - 


Hw, G0 fi" (n)sign(f, M) (n) + 
; 1)= w; 21 
HUM | n, GO si) | CP 
Where step size parameter Hw; (n) is calculated as: - 
Hw 
hy) = Q2) 


ll£& G0 $2 G0l[ 


My is experimentally determined by extensive simulations. 
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Fig. 9: Proposed method — 4, NSS MFxLMAT & VSSLMS 


IV. COMPUTATIONAL COMPLEXITIES 


Any algorithm's implementation in real-time applications significantly depends upon computational complexity. 


Table-1 shows summary of computational complexities for various existing and proposed methods discussed in this paper. 


Table-1. Computational Complexities 


Methods x, +, Y +,- Total 
Eriksson's method 16M+8Lw+4 18M+6Lw-6 34M+14L,,-2 
Akhtar’s method 16M+8L,,+22 18M+6L,,+6 34M+14L,,+28 
Ahmed’s method 24M+14L,,+8 1 28M+14Ly+34 52M+28Ly+115 
Yang’s method 18M+8Lw+4K+40 20M+6L,,+4K+12 38M+14Ly+8K+52 
"— 4L,7+4M?+24M+ ALs^-4M?-18M^ 8L,28M?-42M-4- 

14Lw+4K+34 6Lw+4K+12 20Lw+8K+46 
Proposed method — 1 ALs^-16M--14L 422 AL y?+18M+6Ly+6 8Ly7+34M+20LW+28 
Proposed method — 2 16M+8L,,+26 18M+6L,+10 34M+14L,,+36 
Proposed method — 3 16M+8L,,+30 18M+6L,,+6 34M+14L,,+36 
Proposed method — 4 2L+16M+6L,+36 4L+18M+6L,+4 6L4+34M+12L,+40 
V. SIMULATIONS AND DISCUSSIONS performance of proposed methodologies by comparing 


MATLAB simulations carried out to validate the superior 


them with existing methods as mentioned below: 


>  Erikson's method [11] 
>  Akhtar's method [12] 


This section discusses results of extensive 
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»  Ahmed's method [13] 
> Yang's method [14] 
»  Jabeen's method [16] 
Mean Noise Reduction (MNR;) and Relative 
Modeling Error (AS;;) are two performance metrics used 


for comparison of all algorithms under investigation. The 
MNR; is defined as: - 


(23) 


MNR,(n) = E ts] 


Ag k (n) 
where Ae, (n) is measurement of absolute value of 
disturbance signal and Ag, (n) is absolute value of residual 


error, both measured at k" error microphone, calculated as: 


Ae, Qn) = Ag (n - 1) + A — Alex) (24) 
Ag,() = AAg, (2-1) + (1 - 2)|d Qn) (25) 
On the other hand, the value of A5,; is calculated as given 
in (26) 
lls G0 — su|] 
Ils; Qo 


Using data from [4], original primary P,;(n) and 


AS,;(n) = 201ogio 


(26) 


secondary S,;(n) acoustic paths are taken as FIR filters. 


Impulsive Input with a*1.85 
T T 
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Values of various fixed parameters used for simulations in 
this paper are given in Table-2. Extensive simulations have 
been carried out to determine appropriate values for various 
controlling parameters to attain best results and these values 
are given in Table 3. Modeling filters S, j (1) in all proposed 
methods are adapted through VSSLMS algorithm which has 
stability issue at start when step size is smaller [31]. To 
avoid instability, modelling filters 5, j (0) are initialized by 
offline modelling until the modeling error is dropped to -5 
dB instead of null vector [32]. Cases 1-3 of this section 
discuss performance of proposed algorithms under 
stationary environment with varying impulsive input 
(Error! Reference source not found.). Case — 4 presents 
performance for non-stationary environment: - 


» Case — 1: Moderate impulsive input (a=1.85) 
» Case — 2: High impulsive input (a=1.65) 

» Case — 3: Very high impulsive input (a=1.45) 
» Case — 4: Non-Stationary Environment 


(1 guess you have used subplotting command, don't use it, 
also pics are still very unclear) if you want to plot more than 
one pic and preseny them simultaneously, plot them 
separately, insert in word, label a,c\b,c,d, etc and add 
captions for each of them at end 


Impulsive Input with o=1,65 
T T 


2 
S 
z 
S 


a 
Ss 
S 


E 3 
> E 
a a % 
E A 
j i -500 
1000 
Sample Number x104 Sample Number x10* 
40000 Impulsive Input with mts 40000 — . i reed Impulsive Input — i , 
8000 8000 
3 6000 3 co 
S 4000 3 
H 2000 E iii 
4 Ü dj 2000 
-2000 0 
-4000 — -2000 — " 
1 2 3 4 5 6 7 8 9 10 
Sample Number x10* Sample Number x10* 
Fig. 10: (a) Moderate impulsive input (a=1.85), (b) High impulsive input (a=1.65), 
(c) Very High impulsive input (a=1.65), (d) Blended impulsive input 
Table-2. Various simulations parameters 
MC ANC system with OSPM IN [x,. (7) ] 
Parameter Symbol Value Parameter Symbol Value 
Primary paths tap size L 48 Total samples N 100,000 
Secondary paths tap size M 16 Total realizations Avg 10 
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Control filters tap size Ls 32 Characteristic exponent " 1.85, 1.65, 1.45 
Y 
OSPM filters tap size M 16 Other parameters 0 
ô 
Table 3. Controlling Parameters 
Methods Case — 1 Case - 2 Case — 3 Case- 4 
Uw = 5e-7 
Eriksson's method Ls = le-4 - - - 
A= .99 
Uw = 3e-7 
Akhtar's method Hired - - - 
[ls max = 7.5e-4 
A= .99 
fy = 2.5e-3 
u = 9.5e-3 
Ahmed’s method dcm - - - 
Ymin = .3 
Ymax = .9 
A= .999 
lw = le-7 Uw = le-7 Ly = le-8 Ly = le-8 
Tar einh un = le-4 un = le-4 Ln = le-8 Ln = le-8 
a = .005 a = .0025 a= .0015 a= .0015 
A= .9999 A= .9999 A= .9999 A= .9999 
ó; = le5 ó; = le5 ó; = 1e6 ó; = 1e6 
T ANNE 62 = 4.5e4 02 = 6.5e4 à? = 6.5e5 62 = 6.5e5 
un = le-4 un = le-4 Ln = le-6 Ln = le-6 
A= .99 A= .999 A= .99 A= .99 
ô = 5e4 ô = 5e4 6 — 5e5 6 = 5e5 
Proposed method - 1 Beg Tee M bere Heme de 
Is max = Te-3 s. max = Te-3 Ils max = Te-3 Ls_max = Te-3 
A= .99 A= .99 A= .99 A= .99 
Lw_min = 1e-6 Lw_min = 1e-6 [Iw min = 1e-7 [Iw min = 1e-7 
[Aw max = 1e-4 [Aw max = 1e-4 Hw max = 1e-5 [ls max = 1e-5 
Proposed method - 2 Ils. min = 7.5e-4 Ils min = 7.5e-4 Ils min = 7.5e-6 [ls min = 7.5e-6 
ls max = 7.5e-3 ls max = 7.5e-3 ls max = 7.5e-3 ls max = 7.5e-3 
A= .99 A= .99 A= .99 A= .99 
ee ee ee [Ay = le-7 Ly = le-9 Uw = 5e-13 Uw = 5e-13 
Lw_max = 7e-2 Lw_max = 7e-2 [hw max = 7e-2 Lw_max = Te-2 
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Ils. min = Le-4 Ils min = le-4 Is min = le-6 Ils min = le-6 

A= .999 A= .999 A= .999 A= .999 

Uw = 5e-3 Uw = 5e-3 lUw = 9.5e-4 [Ay = 9.5e-4 
Proposed method - 4 Hunde aa La Bye ces 

s max = Te-3 Is max = 7e-3 Is max = 7.5e-3 [ls max = 7.5e-3 

À-.99 À-.99 A= .99 A= .99 
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Case — 1 : Moderate Impulsive Input (a—1.85). To begin 
with simulation part, all reported algorithms [11, 12, 13, 14, 
16] are subjected to moderated impulsive noise (a=1.85) for 
MC ANC with OSPM. Fig. 11 shows comparison of relative 
modelling error AS,;(n) and mean noise reduction 
(MNR,(n)) of various existing techniques. Fig. 11 (a) 
shows that Yang's method exhibits fastest convergence at 
n=10,000 and achieving lowest value of AS;;(n) = 


—25 dB while MNR curve depicted in Fig. 11 (b) shows 
that Jabeen's method achieved convergence at n=20,000 
with lowest steady state value (0.35 dB) among all the 
existing methods [11, 12, 13, 14, 16]. Since, Yang's and 
Jabeen's methods outperform other exiting methods hence, 
only these two methods are being employed in future 
simulations and compared with algorithms proposed in this 
paper. 


Relative Modeling Error : Case-1 (Impulsive Input with a=1.85) 
o 


iterations x10* 


| Eriksson's Method 


= = -AKkhtar's Method 
— — Shakeel's Method 
Yang's Method 

|—— Faresha's Method | 


A S21 [dB] 
S à 


4822 [dB] 


iterations x10* 


77 Eriksson's Method. 

= = = Akhtar's Method 

= Shakeel's Method 
Yang's Method 


Fareeha's Method 


iterations 


Fig. 11: Simulation Results - Existing Methods, Case-1. (a) AS (b) MNR 
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Although, Jabeen's method [16] exhibited best 
performance in terms of noise reduction among existing 
techniques as shown in Fig. 11 (b) however, this 
exceptional performance was achieved at the cost of higher 
computational complexity due to use of FxRLS-FxRLS 
algorithm in both, control W;(n) as well as modelling 
S (n) filters. To overcome this additional computation 
load while achieving similar performance, we developed 
our first proposed method which is a combination of 
FxRLS-VSSLMS algorithms (Proposed method - 1). This 
proposed method has quite less computational complexity 
(Table-1) as compared to Jabeen's method but in the noise 
reduction, it even surpasses Jabeen's method in terms of 
convergence speed, robustness to impulsive input and lower 
steady state value (Fig. 12). MNR curve depicted in Fig. 12 
(b) shows that Proposed FxRLS-VSSLMS achieved 
convergence at n=8000 with steady state value of 0.27 dB. 


Good results of Proposed FxRLS-VSSLMS 
algorithms motivated us to further develop our second 
proposed algorithm, VSSLMS-VSSLMS (Proposed method 
- 2), which has even lesser complexity (Table-1) than our 
proposed method 1 and still achieves same steady state 
MNR values (0.27 dB) as in the case of Proposed method 1, 


e$ h e 
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however, with reduced convergence speed at n = 16,800 as 
shown in Fig. 12 (b). 


To combine the positive aspects (e.g. low steady 
state error and robustness) in both the proposed methods — 
1 and 2 and cater for their negative aspects (e.g. high 
computational complexity in Proposed FXRLS-VSSLMS 
and slow convergence in Proposed VSSLMS-VSSLMS), a 
third method is proposed. This third proposed method 
combines FxLMAT (as used by Khan et al. against IN [30]) 
with VSSLMS algorithms (Proposed method — 3). This 
method when implied in a MCANC system in the presence 
of IN, although achieved steady state value (0.35 dB) as in 
Jabeen's method, however, its convergence is very slowly, 
achieved at n=50,000. Moreover, Proposed FXLMAT- 
VSSLMS does not show robustness to impulsive nature of 
input noise, (Fig. 12 (b)). Hence, lastly, a modified version 
of FxLMAT is presented in proposed method 4 as 
combination of NSS MFxLMAT-VSSLMS. Fig. 12 (b) 
shows that this proposed NSS MFxLMAT-VSSLMS 
demonstrates best results among all techniques with fastest 
convergence at n=4000 and same steady state value of MNR 
= 0.27 dB with least complexity (Table-1). 


Relative Modeling Error : Case-1 (Impulsive Input with «=1.85) 
at a o [———————————4— 


iterations x 10^ 


iterations x10* 
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MNR : Case-1 (Impulsive Input with c1.85) 
T T T 
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iterations 


Fig. 12: Simulation Results, 


Case — 2: High Impulsive input (a—1.65). For high IN, 


similar results are obtained as in Case-1 of moderate IN. 
Simulation results shown in Fig. 13 (a) confirms that Yang's 
method [14] extended to MC again gives best result in 
reduction of relative modelling error AS(n) = -30 dB while 
proposed algorithms yet again give best MNR curve but 
with less computational complexity (Table-1). It is evident 


AS11 [dB] 


" i 4 " 1 -L " ES n 
0 1 2 3 4 5 6 7 6 9 10 
iterations : ——Àá44n$ 
Yang's Method 
Fareoha's Method 


A821 [dB] 
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L] 1 2 3 4 5 L] 7 B 9 10 
iterations «104 


Case — 3: Very High Impulsive input (g=1.45). In this 


case, when the input noise becomes highly impulsive, 
Yang's method [14] consistently performed best (Fig. 14 
(a)) for reduction of relative modelling error AS(n)- -30 
dB. On the other end, MNR performance of Proposed 
VSSLMS-VSSLMS has deteriorated for this excessive IN 
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Case-1. (a) AS (b) MNR 


from Fig. 13 (b) that proposed methods — 1, 2 & 4 achieve 
lowest steady state MNR value of 0.24 dB. Proposed 
VSSLMS-VSSLMS exhibits slower convergence at n= 
13500, similar to its performance in case — 1. Proposed 
FxRLS-VSSLMS and Proposed NSS MFxLMAT- 
VSSLMS show faster convergence at n= 8000 and n=4000, 
respectively. 


Relative Modeling Error : Case-2 (Impulsive Input with (71.65) 
T [] T 


A812 [dB] 


iterations x10* 


" n 
0 1 2 3 4 5 6 7 8 9 10 
iterations x10* 


(Fig. 14 (b)). However, Proposed FxRLS-VSSLMS and 
Proposed NSS MFxLMAT yet again performed 
unswervingly best for MNR curve with fastest convergence 
at n=8000 and lowest steady state error of 0.22 dB and 0.28 
dB respectively. 
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MNR : Case-2 (Impulsive Input with a=1.65) 
T T T 
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Fig. 13: Simulation Results, Case-2. (a) AS (b) MNR 


Relative Modeling Error : Case-3 (Impulsive Input with o=1,45) 
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Fig. 14: Simulation Results, Case-3. (a) AS (b) MNR 
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Case — 4: Non-Stationary Environment. In our final 
evaluation, proposed methods are subjected to non- 


stationary environment to check for their robustness as 
follows:- 


a. In first instance, all secondary paths are perturbed 
at iteration n=50,000 and very high IN (o=1.45) of 
case 3 is used as input to assess behavior of 
proposed algorithms. 

b. In second trial, a blended impulsive input 


comprises of varying o is used as under Fig. 10 


(d):- 


4 


1.65 
1.45 
1.85 


0 «n < 35,000 
; 35,000 € n « 70,000 
; 70,000 € n < 100,000 
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Simulation results (Fig. 15 & Fig. 16) confirm 
robustness of proposed algorithms under non-stationary 
environment.Fig. 15 (a) shows that only Yang's method, 
Proposed FXRLS-VSSLMS and NSS FxLMAT-VSSLMS 
converge after perturbation is encountered at n=50,000. 
Moreover, Yang's method consistently manifested lowest 
relative modeling error by reaching to -30 dB. Similarly, 
Fig. 15 (b) shows that Proposed FxRLS-VSSLMS and NSS 
MFxLMAT-VSSLMS performed in same manner as in case 
— 3 (a=1.45). In this case, both the proposed algorithms have 
touched lowest steady state error of 0.20 dB and 0.22 dB 
respectively. It is important to note that after encountering 
perturbation at n=50,000, both algorithms took 8000 
iterations again to converge at n=58,000 as in the beginning 
(n=0 to n=8000). 


Relative Modeling Error : Case-4 (Non-Stationary Environment for Impulsive Input with n1.45) 
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Fig. 15: Simulation Results, Case-4 (a — 1.45). (a) AS (b) MNR 
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Relative Modeling Error : Case-4 (Non-Stationary Environment for blended Impulsive Input) 
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Fig. 16: Simulation Results, Case-4 (blended input). (a) AS (b) MNR 


Fig. 16 shows response to blended IN. At 
beginning n=0, high impulsive input (0=1.65) is used. At 
n=35000, it is changed to very high IN (a=1.45) and at 
n=70000, to moderate IN (a=1.85). It is evident from 
simulation results depicted in Fig. 16 (a) that Yang’s 
method continue to perform best for relative modeling error. 
A slight change is observed in convergence pattern at swap 
over points in relative modeling error. Fig. 16 (b) illustrates 
MNR for blended IN. Proposed FxRLS-VSSLMS and NSS 
MFxLMAT-VSSLMS once again reached to lowest steady 
state error of 0.22 dB and 0.25 dB, respectively. Only a 
minor disturbance is noticed at swap over points (n=35000 
and n=70000). 
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VI. CONCLUSION 


In this paper, four new methods are proposed to 
mitigate IN in MC ANC system employing OSPM for the 
first time under stationary as well as non-stationary 
environment. The outcomes of the simulations show that the 
proposed algorithms offer quicker convergence and lowest 
steady state error for MNR than existing approaches with 
comparable modeling accuracy. Proposed NSS 
MFxLMAT-VSSLMS achieved same performance as par 
with high order complex FxRLS algorithm with much less 
computational complexity. There is a room of improvement 
available in modeling accuracy which is a task of future 
work. 
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